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CHAPTER  1 


INTRODUCTION 


This  report  is  the  result  of  CENSEI's  basic  research  work  on  critical 
issues  concerning  system  engineering  and  integration  of  spread  spectrum 

systems,  in  ooncert  with  applications  of  distributed  systems  architectures,  to 

2 

provide  information  distribution  for  C  information  management  systems 
operating  in  a  communications -bounded  environment.  In  a  communications- 
bounded  environment,  subscribers  on  a  given  network  are  transmitting  at 
maximum  subscriber  capacity.  Network  capacity  is  not  the  simple  sum  of 
subscriber  capacities.  Rather,  it  can  be  a  complex  function  of  the  number  of 
subscribers,  their  topological  distribution,  and  their  offered  traffic  load. 
The  RF  resource,  however,  is  limited,  and  saturation,  which  translates  to 
self-jamming,  will  result  as  more  subscribers  are  operating,  or  as  more 
capacity  is  allocated  and  utilized  by  each  subscriber.  The  long  term  goal  is 
to  derive  a  meaningful  portrayal  of  performance  degradation  due  to  self¬ 
jamming  in  spread  spectrun  system  candidates  for  the  Army's  survivable  C3 
network.  Such  a  portrayal  will  then  be  used  by  the  Army  system  engineers  to 
ensure  that  network  capacity  is  adequate  for  given  user  applications. 
Initially,  performance  degradation  was  studied  for  hybrid  direct-sequence/ 
frequency-hop  systems.  In  the  course  of  this  study,  it  became  apparent  that 
performance  degradation  was  also  a  function  of  the  receiver  decision  logic, 
which  is  the  subject  of  this  report. 

The  spread  spectrum  communications  system  which  was  studied  uses  an 
alphabet  of  32  symbols  to  communicate  messages.  A  symbol  is  represented  by  a 


sequence  of  32  chips.  (A  chip  is  a  binary  (+1,-1)  digit  which  takes  on  a 
value  of  either  +1  or  -1.)  A  sequence  of  32  chips  is  cyclically  shifted  to 
generate  the  entire  32-symbol  alphabet. 

The  system  uses  a  Reed -Solomon  code  for  channel  encoding.  (The  reader  is 
referred  to  Macwilliams  &  Sloane  [1]  for  a  discussion  on  the  Reed-Solomon 
codes.)  An  R-S  (31,15)  code  is  used  to  encode  15  information  symbols  into  31 
channel  symbols,  so  that  16  symbols  are  used  for  error  correction  and  error 
detection  purposes.  The  code  allows  for  forward-error  correction  if 

2e  +  E  <  17  (1-1) 

where  e  =  the  number  of  symbol  errors  and  E  «  the  number  of  symbol  erasures 
within  the  31-symbol  code  word.  If 

2e  +  E  2  17  (1-2) 

then  a  decoding  error  will  occur. 

The  probability  of  a  word  error  (Pw),  or  the  expected  relative  frequency 
of  a  word  decoding  error,  vrtiere  a  word  is  represented  by  an  R-5  (31,15)  code, 
can  be  formulated  using  the  probability  of  a  symbol  error  (Pg)  and  the 
probability  of  a  symbol  erasure  (P£)  as  independent  variables.  But  PE  and  P@ 
are  not  independent  variables:  they  are  dependent  upon  the  signal-to-noise 
ratio  (SNR),  the  modulation  technique  used,  and  the  threshold  setting  for  the 
decision  rule  used.  Nevertheless,  constant  P.,  curves  as  a  function  of  P_  and 
P@  can  be  plotted  on  log  x  log  graph  paper.  The  set  of  Pw  curves  becomes  a 
powerful  tool  vrtien  used  in  conjuction  with  the  P_  vs.  P  plot  with  constant  SNR 


curves,  because  by  superimposing  the  curves  we  can  relate  SNR  to  P^. 

Decision  rules  for  detection  of  spread  spectrum  signals  were 
investigated.  Decision  rules  were  designed  using  simple  and  ratio 
thresholding  techniques.  For  each  decision  rule,  the  performance  of  the 
detection  system  in  terms  of  P_  and  P  was  analyzed.  The  obtained  data  were 
used  to  plot  constant  SNR  curves  and  constant  threshold  curves.  With  the  use 
of  the  P^  curves  the  optimum  threshold  setting  was  found  for  a  given  SNR. 

Orthogonal  signalling  (FSK)  and  antipodal  signalling  (PSK)  were 
investigated  to  determine  their  effects  on  the  performance  of  the  detection 
system. 

Analog  and  digital  correlations  were  considered.  The  performance  curves 
of  analog  correlation  and  digital  correlation  were  compared. 

Monte  Carlo  simulation  of  the  digital  detection  process  using  the 
threshold  decision  rule  was  used  to  produce  the  performance  data  which  can  be 
compared  with  the  results  of  the  analysis. 


CHAPTER  2 


DETECTION  THEORY  APPLIED  TO  32-ARY  ORTHOGONAL  SIGNALS 

For  radio  communications,  a  transmitted  signal  becomes  corrupted  with  an 
additive  gauss! an  white  noise  when  it  reaches  its  destination.  Because  a 
signal  attenuates  during  its  propagation,  signal  degradation  due  to  the 
presence  of  the  noise  in  the  signal  becomes  a  problem  especially  for  long 
distance  transmissions.  Because  the  white  gauss i an  noise  is  a  random  process 
the  likelihood  of  a  detection  error  occurring  is  described  in  terms  of  the 
probability  of  an  (detection)  error. 

A  received  signal  is  represented  by 

y(t)  -  x(t)  +  n^t)  (2-1) 

where  x(t)  is  the  transmitted  signal  and  nw(t)  is  the  gaussian  white  noise. 
The  noise  nw(t)  is  assumed  to  have  the  white  spectral  density 

G(f)  »Nq/2.  (2-2) 

The  receiver  operates  cm  the  received  signal  y(t)  to  determine  which  of  the 
symbols  i*l,2,...,H,  were  sent. 

If  is  a  symbol  belonging  to  the  set 

U  *  { Uj ,  •••»  )^|}>  (2—3) 

then  yj  is  an  M-ary  symbol  that  is  equivalent  to  log2(M)  bits  of  information. 


The  spread  spectrum  communication  system  in  our  study  uses  32-ary  symbols, 
each  symbol  representing  5  bits  of  information;  32  chips  are  used  to  generate 
a  signal  Xj  (t)  representing  pij  (for  i=l,2,...,32). 

32 

Thirty- two  chips  are  capable  of  generating  2  different  signals; 

32 

however,  among  the  2  prospective  signals,  only  32  which  are  orthogonal  to 
each  other  are  chosen.  There  are  various  ways  of  generating  32-ary  orthogonal 
signals,  and  we  will  explore  two  of  the  methods  in  Chapter  6. 

Supjpose  there  exists  a  set  of  32  orthogonal  signals  described  by 

X  ■  (x^t),  x2(t),  ...,  x32(t)}.  (2-4) 

A  set  of  vectors 

V  =  {^(t),  4>2(t) ,  ...,  «t>32 (t) }  (2-5) 

forms  the  basis  of  the  vector  space  s32  of  dimension  32  which  is  spanned  by 
the  linearly  independent  vectors  x1(t),  x2(t),  ...,  x32(t).  Each  signal  Xj(t) 
is  a  linear  combination  of  <t>j,  <P2,  ...,  $32  such  that 

xi(t)  "  xil*l  +  xi2*2  +  —  +  xi  32^*32  (2_6) 

where  x^,  Xj2»  ...,  x^32  are  scalars.  The  orthogonality  of  the  signals 
implies  that 


=  Es 

1=1,2,  . . . , 32 

(2-7a) 

<Xj,Xj>  =  0 

if). 

(2-7b) 
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2.1  Signal  Space  Formulation 


The  present  discussion  of  signal  space  formulation  is  limited  to  the  case 
of  having  32-ary  orthogonal  signals.  The  reader  is  referred  to  Carlson  [2J 
for  a  generalized  treatment  of  signal  space  formulation. 

Suppose  there  are  32  orthogonal  signals  Xj(t),  x^it),  ...,  x32(t).  Then 
there  exists  a  signal  space  c32  containing  all  of  the  above  signals;  C32  *s 
spanned  by  an  orthonormal  basis  {^(t),  4>2(t),  ...,  ♦32(t)},  so  that 

*ilt)  -fj  xiA<t>  <2-« 


An  additive  zero-mean  gaussian  white  noise  present  in  the  transmission 
medium  may  not  be  fully  contained  in  c32'  in  which  case  it  cannot  be  fully 
described  by  a  linear  combination  of  the  orthonormal  basis  vectors  ♦jCt), 
<j>2(t),  ...»  <j>32(t)  spanning  the  signal  space  532.  But  the  noise  nw(t)  can  be 
described  as  a  sum  of  two  terms,  one  term  representing  the  component  vector 
which  does  belong  to  C32  and  the  other  term  representing  the  component  vector 
which  does  not  belong  to  ?32.  The  noise  can  be  represented  as 

n  (t)  *  n(t)  +  n  (t)  (2-10) 

w  e 

where  n(t)  is  the  projection  of  the  vector  n,,(t)  on  C,-  and  n  (t)  is  the 


irrelevant  noise 


The  relevant  noise  n(t)  is  a  linear  combination  of  the  basis  vectors  and 
is  described  by 


32 

n(t)  =  £  i  nkVt)  <2-11) 

where 

n^t)  -  <n,4^>  =  <VV*  (2-12) 

Whereas,  for  the  irrelevant  noise, 

^e'V  *  0  k«l,2,...32  (2-13) 

because  ne(t)  does  not  belong  to  the  signal  space  c32» 

The  detector  for  the  system  in  our  study  has  a  bank  of  32  correlators, 
received  signal  is  correlated  with  32  seperate  locally  generated  signals  so 
that 

<y,x.>  =  r. 

<y'x2>  ■  r2 
•  • 

•  • 

<y,x32>  -  r32  (2-14) 

where  rk  is  the  output  of  the  k’th  correlator. 

Note  that 

y(t)  «  Xj(t)  +  nw(t)  *  Xj (t)  +  n(t)  +  ne(t)  (2-15) 


Therefore,  substituting  for  y(t). 


<y,x.>  *  <x.+n+n  ,x.> 

J  1  ®  J 

<y#x.>  *  <x,,x.>  +  <n,x.>  4-  <n  ,x.>. 

J  *  D  J  ®  D 


<vy 


<y#Xj>  *  +  <n,x^> 

<y,Xj>  *=  <xi+n,Xj>. 


Defining  z(t)  as  the  projection  vector  of  y(t)  on  C32» 
z(t)  «  Xjit)  +  n(t) 


<y»*j>  *  <z#Xj>. 


(2-16a) 
(2-1 6b) 


(2-17) 


(2-18a) 
(2-1 8b) 


(2-19) 


(2-20) 


Assuming  that  only  hard  decisions  are  made,  the  symbol  corresponding 
to  the  signal  x^(t)  to  which  z(t)  lies  closest  in  the  signal  space  C32  *s 
chosen.  For  z(t)  to  lie  closest  to  x^(t)  the  following  condition  must  be 
satisfied: 


|z-Xj|  -e  Jz-Xjl 


all  i^j. 


(2-21) 
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If  Xj(t)  is  the  transmitted  signal,  then 

z(t)  -  Xj (t)  *  n(t) 


and 


(2-22) 


(2-*j||  -  W-  <J-23> 

where 

j|njj  -  <n,n>1/2.  (2-24) 

The  inner  product  of  the  relevant  noise  with  itself  is  the  sum  of  the 
magnitude  squared  of  the  component  vectors  as  follows: 

<n,n>  »  <n14>1+n24)2+  •••  +n32<*>32'nl<*>l4n2<*>2+  ***  +n32<**32>  (2-25a) 

<n,n>  «  n2  +  n2  +  ...  +  n22  (2-25b) 


where  each  n^  has  the  following  statistical  variance: 

,2 


V  -  "ni  -  E<niU 

1  2  — 5  _  2 
ani  =  ni  "  ni  * 


But 


and 


ni  *  <nw'V 


(2-26a) 
(2-2  6b) 


(2-27) 


(2-28) 
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because  nw(t)  is  a  zero-mean  gauss i an  white  noise.  Therefore, 


(2-29) 


ni 


0. 


Using  (2-29),  (2-26b)  is  simplified  as 


(2-30) 


-  *  2 
The  term  n^  ,  which  can  be  written  as  E[nj  ],  is  analyzed  by  first 

deriving  the  relationship  for  Etn^].  Etn^n^]  is  defined  as 

Etnjn.]  «  e[j[^  nw(t)^(t)dt£o  Vmm(1)dl]  .  (2-31) 

But  since  only  nw(t)  and  nw(l)  are  random  variables  Etn^]  can  be  rewritten 


Efn^J  «  B(nw(t)nwa)J#i(t)^(l)dtdl.  (2-32) 

Efn  (t)n  (1)]  can  be  represented  as  an  inner  product  of  n  (t)  and  n  (1)  as 

W  W  W  w 

Elnw(t)hw(l)l  *  <nw(t),nw(l)>  (2-33) 

E[nw(t)nw(l)l  *  <nw(t),nw(t-(t-l))>.  (2-34) 

i 

Noting  that 


R  (t-1)  «  <nft),n  <t(t-l))> 
nw 

a  new  relationship  for  E[nw(t)nw(l)l  is  found  as  follows: 

Efn  (t)n(l)]  *  R  (t-1). 
w  w  nw 


(2-35) 


(2-36) 
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The  autocorrelation  of  white  noise  is 


(2-44) 


P_(n)  *  Pn  ^ni,Pn  <ny^> 
n  n  j  1  rij  i 


P_  (n,_) 

n32  32 


to  get 

Pn(n)  *  ( 7JNq  ) _1 6exp  [ -  |n || 2 A* Q  J 

(2-45) 

where 

|jnjj2  *  n^2  +  r^2  +  ...  +  n22^  * 

(2-46) 

I 

I 


2.2  Maximum  Likelihood  Decision  Rule 


Given  a  projection  vector  z(t),  it  is  possible  to  choose  the  signal  most 
likely  to  have  been  transmitted.  We  impose  the  condition 

P(Hj|z)  >  PfWjJz)  all  i^j  (2-47) 

which,  if  satisfied,  allows  us  to  decide  that  y..  is  the  most  likely  symbol  to 
have  been  transmitted. 

Using  the  Bayes'  rule 

Piyjz)  *  p(yi)P2(z|yi)/P2(z)  (2-48) 

(2-47)  is  restated  as 

P(Vj>P2<*|l»j)  >  P(yi)P2(z|Wi)  all  i*j.  (2-49) 

Assuming  that  every  symbol  has  an  equal  probability  of  being  transmitted, 
(2-49)  is  simplified  as 

Pg(s|»j)  >  Pz(z|V  all  i#j.  (2-50) 

Since  z(t)  given  Uj  implies 

n(t)  =  z(t)  -  Xj(t)  (2-51) 

where  is  represented  by  Xj(t),  pz(z(yj)  can  be  written  as 

Pz(z|wi)  -  Pp(z-xi).  (2-52) 

Using  (2-45)  and  (2-52)  we  get 


I 


I 


VEIV  =  («0)~16exp[-^-x.||2/N0I 
and  using  (2-50)  and  (2-53)  we  get 

expC-Jz-XjH2/^]  >  exp[-gz-x.g2/N0] 

Which  is  simplified  as 

I2-*}#2  <  <z~xill 2  aU 

Where 

Uz-xJ2  =  <z-x..z-x.> 

<z-xi,z-x.>  =  Uz||?  -  2<z,x.>  -  jJx.JI2. 

The  energies  for  the  signals  z(t)  and  x^(t)  are  defined  as 
Ez  -  Rz(0>  -  |z||2 

E«i  *  \<0)  -  hil2- 

Using  the  above  definitions,  (2-57)  is  restated  as 

<z-x.,z-x.>  =  E  -  2<z,x.>  +  E 
vi  z  i  x^ 

Which  is  equivalent  *  ~ 

<z-x.,z-x.>  =  E  -  2<y,x.>  +  E 

l  1  Z  l 

because 

<z,xL>  =  <y-ne#x.>  =  <y,x.>. 


11 


[rxr^/, 


Equations  (2-55),  (2-56),  and  (2-60)  are  used  to  get 


<2, x j>  -  (l/2)E..  >  <z,x.>  -  (l/2)E.  all  ifcj.  (2-62) 


But  we  have  a  signal  set  vAxse  32  orthogonal  signals  have  the  property 


Ej  -  Ej  =  ...  *  E32. 


When  (2-63)  is  considered,  (2-62)  sinplifies  to 


<z,x.>  >  <z,x.> 

3  i 


(2-63) 


all  i£j. 


(2-64) 


lb  choose  the  synixol  most  likely  to  have  been  sent,  the  modified  decision 


function  is  defined  as 


Di'  -  <z,xt>  «  <y,xA> 


(2-65) 


and  vij  is  chosen  such  that 


D.'  >  D. ' 
3  i 


all  i£j. 


(2-66) 


2.3  Decision  Rules  Utilizing  the  Erasure  Symbols 

The  detector  in  our  study  is  capable  o£  declaring  symbol  erasures.  A 
symbol  erasure  is  declared  when  the  received  signal  y(t)  has  unacceptably 
degraded  in  signal  quality  due  to  the  noise  present  in  the  transmission 
medium.  Expressed  in  another  way,  a  symbol  erasure  is  declared  when  the 
probability  of  a  symbol  error  becomes  unacceptably  high. 

For  the  case  of  having  a  detector  which  is  capable  of  declaring  erasures, 
the  function  of  the  decision  rule  is  not  to  minimise  the  P ,  but  to  find  the 
optimum  ratio  of  P  to  P  .  Unfortunately,  the  optimum  ratio  is  not  fixed,  it 
is  dependent  upon  SNR,  the  modulation  technique  used,  and  the  decision  rule 
used. 

The  materials  subsequent  to  Chapter  3  will  analyze  the  optimization  of 
detection  performance  using  different  options  in  the  decision  rule  used, 
modulation  technique  used,  and  the  correlation  method  used.  But  we  will  first 
investigate  P„  as  a  function  of  P_  and  P  . 

W  u  6 


CHAPTER  3 


PROBABILITY  OF  A  WORD  ERROR 


The  discussion  in  this  chapter  is  based  on  a  report  fcy  Rosenstark  and 


Frank  [3]. 


In  this  chapter,  the  term  word  error  refers  to  the  decoding  error  of  a 
Reed-Solanon  (31,15)  code.  Thirty-one  detected  symbols,  each  of  which  nay  be 
a  correct  symbol,  a  wrong  symbol,  or  an  erasure,  are  processed  by  the  decoder 
to  generate  a  vord  consisting  of  15  information  symbols.  The  decoded  word, 
however,  can  be  in  error  when  there  exists  enough  symbol  errors  and  erasures. 
This  chapter  studies  the  relationships  among  P  ,  P  ,  and  P  . 

£j  G  W 

The  system  in  our  investigation  uses  the  Reed-Solanon  (R-S)  coding.  The 
R-S  (31,15)  decoder  is  capable  of  error-free  decoding  if 


2e  +  E  <  17 


where 


e  =  the  number  of  symbol  errors, 

E  =  the  number  of  symbol  erasures. 


A  word  (code)  error  occurs  if 


2e  +  E  =  S 


where 


S  2  17. 


(3-1) 


(3-2a) 
( 3-2b) 


(3-3a) 


::-3b) 


l_i  : 


Table  3.1.  Some  Ccntoinations  of  Errors  and  Erasures 
Which  Cause  Code  Word  Errors 


2e+E=17 

2e+E=18 

2e+E=19 

2e+E=20 

e 

E 

e 

E 

e 

E 

e 

E 

0 

17 

0 

18 

0 

19 

0 

20 

1 

15 

1 

16 

1 

17 

1 

18 

2 

13 

2 

14 

2 

15 

2 

16 

3 

11 

3 

12 

3 

13 

3 

14 

4 

9 

4 

10 

4 

11 

4 

12 

5 

7 

5 

8 

5 

9 

5 

10 

6 

5 

6 

6 

6 

7 

6 

8 

7 

3 

7 

4 

7 

5 

7 

6 

8 

1 

8 

2 

8 

3 

8 

4 

9 

0 

9 

1 

9 

2 

10 

0 

Defining  e  .  as  the  minimum  nurrber  of  syntool  errors  required  for  a  word 
min 

error  to  occur,  we  see  from  table  3.1  that 


for  at  least  S  equal  to  17,  18,  19,  .and  20.  We  will  define 
maximum  number  of  symbol  errors  for  a  given  value  of  S.  An 
table  3.1  shows  that 


(3-4) 

e  as  the 
max 

inspection  of 


e 

max 


<  S/2. 


(3- 5a) 


It  is  also  apparent  that 


where  jxj,  the  floor  of  x,  is  defined  as  the  largest  integer  less  than  or 
equal  to  x  (the  definition  appears  in  Iverson  [4]). 

When  a  signal  is  detected,  a  symbol  decision  is  made.  The  decision  can 
be  correct,  be  in  error,  or  be  declared  an  erasure.  The  sum  of  the  number  of 
errors  and  the  nunber  of  erasures  in  an  R-S  (31,15)  word  cannot  exceed  the 
total  number  of  symbols  in  the  word  (there  are  31  symbols  in  an  R-S  (31,15) 
code  word),  or 


e  +  E  <  31 


which  can  be  restated  as 


E  3  31  -  e. 


Using  (3-3a),  the  above  becomes 


E  =  S  -  2e. 


Using  (3-7)  and  (3-8)  we  arrive  at  the  relationship 


e>  S  -  31. 


In  the  above,  e  must  be  a  non-negative  integer,  therefore,  em^n  is  defined  for 
two  different  cases  as  follows: 


S  -  31 


if  S  >  31 
if  S  £  31. 


(3-10) 


Referring  to  (3-6),  the  maximum  possible  value  of  e  is  31.  By  substituting  31 


for  e _ in  (3— 5b)  we  get 

max 

max{S)  =  2(31)  =  62. 

For  a  word  error  to  occur,  S  must  be  an  integer  that 
relationship: 

17<S<  62. 


The  probability  of  a  word  error  is  defined  as 


trtiere  P(S)  is 


w 


62 

s 


P(S) 


17 


P(S) 


max 

P(sne*k). 


"min 


Since  e  .  in  (3-10)  exists  for  two  different  cases, 
mm 

different  cases  as  follows: 


(  S/2 

E  P(S  n  e=k) 
k=S-31 


P(S)  =  \ 


S/2 

£  P(S  fl  e=k) 
k=0 


if  S=*31 


if  SS31 


where  substitutions  are  made  for  e  ,  and  e  ,  using 

max  min 


(3-11) 

takes  on  the  following 

(3-12) 

(3-13) 

(3-14) 

P(S)  also  exists  for  two 

(3-15) 

( 3-5b)  and  (3-10). 


P(SOe=k)  is  defined  as 


(3-16) 


i 


where 


Pg  =  probability  of  a  symbol  error. 

(3-1 7a) 

; 

Pg  =  probability  of  a  symbol  erasure. 

(3-1 7b) 

Pc  =  probability  of  a  correct  symbol. 

(3-17c) 

Noting  that 

r  < 

/31-S+k \  , 

U-S-hJ  * 

(3-18) 

.  4 

(3-16)  is  simplified  as 

«■«-*>•  ft)  (25)  -SW*. 

(3-19) 

m*  ■  »Hi 

But 

:  : 

/31\  /31-k\  _  31!  ( 31— k.) ! 

\  k)  \S-2k /  (31— k) lk!  (31-S+k) ! (S-2k) ! 

(3-2 0a) 

w  *■ 

- 1 

1 

/31\  /31-kN  _  311 

\  k)  \S-2k/  k! (S-2k) i (31-S+k) . 

(3-20b) 

and 

Pc  =  1  -  (Pe+PE). 

(3-21) 

Using  (3-19),  (3-20b),  and  (3-21)  we  get 

— - . 

P(S(1  e=k)  =  3l!P^-2k(l-(Pe+PF;)]31-B+k 

k! (S-2k) ! (31-S+k)! 

(3-22) 

21 


Using  (3-13)  and  (3-15)  we  write  P^  as 

31  S/2  62  S/2 

P  «  Z  Z  P(Sfle*k)  +  Z  Z  P(Sne«k)  (3-23) 

w  S®17  k»0  S=32  k=S-31 

where  P(Sfle“k)  is  defined  by  (3-22). 

p  as  a  function  of  P„  and  P  can  be  numerically  analyzed  using  (3-22) 
w  E  e 

and  (3-23) .  The  obtained  data  can  be  plotted  in  the  form  PE  vs.  Pft  for 

constant  P  .  This  was  done  for  the  values 
w 

Pw  -  Iff2,  10"3,  ...»  10~10  (3-24) 

in  fig.  3.1.  The  Pw  curves  can  be  used  as  a  tool  to  compare  the 
performances  of  various  decision  rules  under  various  detector  operating 


conditions 


CHAPTER  4 


ANALOG  SIGNAL  DETECTION 


Figure  4.1.  Diagram  of  detector  and  decoder. 


Figure  4.2.  32  correlator  detection  system. 


Figure  4.1  illustrates  the  detection  and  decoding  processes  considered  in 
this  chapter.  The  noise  corrupted  signal  y(t)  is  processed  by  the  bank  of  32 
correlators  illustrated  in  fig.  4.2.  Each  correlator  produces  a  correlator 
output  r^  vdiich  is  sent  to  the  decision  logic  device.  The  device  takes  the  32 
correlator  outputs  and  decides  upon  a  symbol  belonging  to  the  set 
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u  -  {V0,  yr  .... 


H31,  erasure} 


(4-1) 


Hie  decoder  receives  the  symbols  until  the  31  symbols  making  up  a  code 
word  are  received.  Forward-error  correction  is  performed  on  the  R-S  (31,15) 
code  word,  and  the  15  information  symbols  are  outputted. 

TWo  different  decision  rules  were  developed  and  analyzed,  the  decision 
rules  studied  are  the  threshold  rule  and  the  ratio  rule. 


4.1  Threshold  Rule 

One  variation  of  the  threshold  rule  works  as  follows:  If  only  one 
correlator  output  equals  or  exceeds  the  threshold,  then  declare  the  symbol 
corresponding  to  the  largest  correlator  output  as  the  received  symbol.  If  the 
above  condition  is  not  met,  then  declare  an  erasure. 

This  decision  rule  has  a  major  flaw  because  the  condition  that  only  one 
correlator  output  must  equal  or  exceed  the  threshold  for  a  symbol  decision  to 
occur  causes  the  likelihood  of  an  erasure  occurring  to  substantially  exceed 
the  likelihood  of  a  symbol  error  occurring  for  the  SNR  ranges  of  our  interest. 
After  reviewing  the  results  of  an  initial  numerical  analysis,  the  analysis  of 
this  decision  rule  was  discontinued. 

The  second  version  of  the  threshold  rule,  which  later  proved  to  yield  a 
proper  balance  of  symbol  erasures  and  symbol  errors,  works  as  follows:  If  one 
or  more  correlator  outputs  equal  or  exceed  the  threshold,  then  declare  the 


symbol  corresponding  to  the  largest  correlator  output  as  the  received  symbol. 
If  either  the  largest  correlator  output  does  not  exceed  the  threshold,  or  if 
we  have  a  draw,  then  declare  an  erasure. 

Note  that  the  threshold  setting  is  variable.  Its  variability  serves  to 
control  the  ratio  of  P^  to  P@.  If  the  incoming  signal  is  weak,  or  if  the  SNR 
is  low,  then  the  threshold  should  be  set  relatively  low  to  prevent  an 
excessive  occurrence  of  symbol  erasures.  But  if  the  incoming  signal  is 
strong,  or  the  SNR  is  high,  then  the  threshold  should  be  set  relatively  high 
so  that  the  maximum  correlation  values  which  are  relatively  low  can  be 
declared  as  erasures.  Under  the  condition  of  a  changing  SNR,  the  threshold 
setting  must  also  change  to  achieve  an  optimun  word  error  performance.  We 
will  later  see  that  this  is  not  the  case  for  the  ratio  rule.  The  distinct 
advantage  of  the  ratio  rule  over  the  threshold  rule  is  that  its  optimum  ratio 
setting  stays  relatively  fixed  for  the  range  of  SNR  of  our  interest. 

The  performance  for  the  detection  system  using  the  threshold  rule  for  its 
decision  logic  device  is  analyzed  for  orthogonal  signalling  and  antipodal 
signalling.  For  an  informative  discussion  on  the  signalling  techniques  the 
reader  is  referred  to  Pasupathy  [5]. 


4.1.1  Analysis  of  the  Threshold  Rule  Performance  for  Orthogonal  Signalling 

Without  going  in  depth  into  the  correlation  analysis  we  will  state  that 
when  an  FSK  modulated  chip  signal  is  correlated  with  a  reference  FSK  modulated 
chip  signal,  we  get 

'  E  if  the  chips  have  the  same  frequency 

r  *  < 

c  0  if  the  chips  have  different  frequencies  (4-2) 

where  rc  is  the  chip  correlator  output  variable  and  Ec  is  the  energy  of  the 
FSK  modulated  chip  signal. 

Suppose  a  noiseless  signal  x^  (t)  is  received.  Because  the  signal 
consists  of  32  modulated  chips,  the  i'th  correlator  peak  output  will  be 

r.  *  32E„  «  E  (4-3a) 

A  C  S 

where  r.  is  the  correlator  output  and  E  is  the  energy  of  the  signal  x. (t). 

X  si 

But  if  Xj (t)  is  correlated  by  a  j’th  correlator  (where  i#j),  then  the 
correlator  output  will  be 


•r3  =  16Ec  =  (1/2)  Es.  (4-3b) 

A  code  word  c.  consisting  of  32  (bipolar)  chips  becomes  a  signal  (t) 
when  modulated.  When  c^  is  correlated  with  c^  (where  i^j),  16  chip  positions 
will  have  the  chips  of  same  polarity  and  the  remaining  lfi  chip  positions  will 
have  the  chips  of  opposite  polarity  (the  reader  can  verify  this  by  using  the 
modified  maximal-length  sequences  in  figure  6.3  or  the  Welti  codes  in  figure 
6.5).  When  c^  is  correlated  with  itself,  all  32  chip  positions  will  have  the 


chips  of  same  polarity.  Using  (4—2)  it  can  be  seen  that  (4 -3a)  and  (4-3b)  are 
true. 

Given  that  every  signal  x^  (t)  (for  i=0,l,...,31)  has  an  equal  probability 
of  being  transmitted,  we  will  assume  throughout  this  report  that  x0(t)  is 
transmitted  without  any  loss  of  general ity- 

Zero-mean  white  gaussian  noise  present  in  the  transmission  medium  is  also 
assumed  throughout  the  paper.  For  orthogonal  signalling,  r0  is  gaussianly 
distributed  around  the  mean  energy  value  of  Eg  with  variance  E^^/2.  The 
other  correlator  outputs  are  similarly  gaussianly  distributed,  but  the 
distributions  have  the  mean  value  of  (1/2) Eg. 

To  get  an  erasure  we  must  have 

ro  <  tEs  I4-4"' 

r,  <  tE  all  <4— 4b) 

where  t  is  the  threshold  variable  and 

0  <  t  2=  1.  (4-5) 

The  probability  of  an  erasure  is  formulated  as 

PE  =  Qt  (l-t)Es/o]  {l-Q[  (t-l/2)Eg/CF)  }31  (4-6) 

where  Q[x]  is  the  probability  that  a  gaussianly  distributed  variable  X  exceeds 
the  value  x  (refer  to  Abramowitz  6  Stegun  [6]). 


A  correct  decision  takes  place  when 


and 


r0  =  y  (4-7a) 

r.  <  y  i*l,2,...,31  (4-7b) 


y  2  tEg.  (4-8) 

The  probability  of  a  correct  decision  is  formulated  as 

<•«  -  jar?  /tLe'  <y‘E31  2/2°2,1-q[  <«-*> 

Three  types  of  a  decision  can  be  made  by  a  decision  logic  device:  a  right 
decision,  a  wrong  decision,  and  a  decision  to  declare  an  erasure.  Regardless 
of  the  probability  of  each  event,  the  probability  of  any  of  the  three  events 
occurring  is  1,  or 

1  ■  P.  ♦  P_  +  (4-10a) 

6  E*  C 

When  Pp  and  P  are  known,  the  probability  of  an  error  is  determined  by  using 

u  C 

the  relationship  in  (4-10a)  as  follows: 

P  *  1  -  P  -  P  .  (4-1 0b) 

6  EC 

Substituting  for  P_  and  P  in  (4-10b)  using  (4-6)  and  (4-9)  we  get 

Ei  C 

Pe  *  l-Qt(l-t)Es/a]{l-Q[(t-l/2)Es/an31 

-y3Fo  /tEge"(y"Es)  /2a  t1-Q[(y-Fs/2)/o]31}dy. 


(4-11) 


The  terms  in  (4-11)  can  be  rearranged  to  yield 


Pe  -  Q[(l-t)Es/a](l-{l-Q[(t-l/2)Es/a]}31] 

"JETo  /tE  e”(y~Es)  /2 °2 1 1- ( 1-Q  f  (yKs/2 )  /a] } 31  ] dy .  (4-12) 

s 

For  very  small  values  of  Q  which  correspond  to  small  values  of  Pg  we  can  use 
the  two-term  binomial  expansion 

(I4x)k  =  1  +  kx  (4-13) 


in  (4-12)  to  obtain 

Pe  «  31Q[(l-t)Es/o]Q[(t-l/2)Es/ol 

+JBTJ  /tE  .•<y’Es,2/2°2Q[  (y-E^2)/o)dy.  (4-U) 

s 

Noting  that 

a  =VESV2  (4-15) 

and 

Es/0  =V2Es/Nq  (4-16) 

(4-6),  (4-12),  and  (4-14)  are  numerically  computed.  With  the  obtained  data, 
PE  vs.  Pg  for  constant  values  of  Es/^0  (SNR)  were  plotted  in  figure  4.3.  The 
loci  of  the  constant  t  are  also  shown.  Superimposing  the  Pw  curves  from 
figure  3.1  on  the  above  curves  enables  us  to  predict  the  detection  system 
performance. 


An  examination  of  figure  4.3  reveals  that  we  should  set  the  threshold  at 
approximately  t*0.8  to  minimize  the  word  error  probability.  Furthermore,  we 
observe  that  the  change  in  word  error  probability  is  minimal  as  long  as  0.7  < 
t  <  0.82,  so  the  threshold  can  change  an  appreciable  amount  before  large 
changes  take  place  in  system  performance. 


Figure  4.3.  Threshold  rule  orthogonal  signalling 
performance  curves. 


lal  Signalli 


Without  going  in  depth  into  the  correlation  analysis  we  will  state  that 
when  a  PSK  modulated  chip  signal  is  correlated  with  a  reference  PSK  modulated 
chip  signal,  we  get 


Ec  if  the  chips  have  the  same  polarity 

-E  if  the  chips  have  the  opposite  polarity 


(4-17) 


where  r  is  the  chip  correlator  output  variable  and  E  is  the  energy  of  the 
c  c 

PSK  modulated  chip  signal. 

Suppose  a  noiseless  signal  (t)  is  received.  Because  the  signal 
consists  of  32  modulated  chips,  the  i'th  correlator  peak  output  will  be 


r.  =  32E  =  E 

i  c  s 


(4-18a) 


where  r ^  is  the  symbol  correlator  output  variable  and  Eg  is  the  energy  of  the 
signal  Xj(t).  But  if  the  j'th  correlator  correlates  the  signal  x.(t)  (where 
i#j),  the  correlator  output  will  be 


(4— 18b) 


A  code  word  Cj  consisting  of  32  (bipolar)  chips  becomes  a  signal  x. (t) 
when  modulated.  When  c^  is  correlated  with  c^  (where  i£j),  16  chip  positions 
will  have  chips  of  same  polarity  and  the  remaining  16  chip  positions  will  have 
chips  of  opposite  polarity  (the  reader  can  verify  this  by  using  the  modified 
maximal-length  sequences  in  figure  6.3  or  the  Welti  codes  in  figure  6.5). 
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When  is  correlated  with  itself,  all  32  chip  positions  will  have  chips  of 
same  polarity.  Using  (4-17),  it  can  be  seen  that  (4-18a)  and  (4-18b)  are 
true. 

The  probability  analysis  for  antipodal  signalling  is  essentially  the  same 
as  the  analysis  for  orthogonal  signalling,  but  the  statistical  mean  value  for 
Tj  (for  i*l,2,...,31)  is  zero  for  antipodal  signalling.  Without  retracing  the 
derivations  we  will  simply  state  the  necessary  relationships  as  follows: 

PE  =  Q[ (l-t)Es/ol {l-Q[tEg/al }31  (4-19) 

Pc  *  JST^/J  e~ ty-y  2/2°2 (l-Q (y/o! ) 31dy  (4-20) 

s 

Pe  «  Qt  (l-t)Eg/a]  [1-U-QftEg/c] }31] 

r  as  2  2 

JtE  e~(y’Es)  /2°  tl-{l“<?[y/o]}31]dy.  (4-21) 

Using  the  definition  in  (4-16),  (4-19)  and  (4-21)  can  be  numerically 
computed.  By  plotting  the  curves  as  was  done  for  orthogonal  signalling  we  are 
able  to  predict  the  probability  of  a  word  error. 

The  curves  derived  from  the  numerical  analysis  are  shown  in  fig.  4.4.  We 
see  that  a  threshold  setting  of  t=0.58  yields  low  word  error  probabilities. 

If  t  is  varied  over  the  range  0.5  <  t  <  0.62  the  effect  on  word  error 
probability  will  be  minimal. 


Figure  4.4.  Threshold  rule  antipodal  signalling 
performance  curves. 


4.2  Ratio  Rule 


The  threshold  rule  was  a  good  decision  rule  to  implement  because  the 
variable  threshold  setting  allowed  for  the  control  over  the  ratio  of  the 
probability  of  an  erasure  to  the  probability  of  an  error.  But  we  noted  that 
the  optimum  threshold  setting  varies  with  the  SNR  of  the  received  signal. 

We  pursued  our  investigation  to  find  a  decision  rule  whose  optimum 
threshold  setting  stays  fixed  when  the  SNR  of  the  received  signal  varies.  The 
decision  rule  that  satisfies  the  requirement  is  the  ratio  rule.  Hie  ratio 
rule  works  as  follows:  Take  the  ratio  of  the  largest  correlator  output  to  the 
next  largest  output.  If  the  ratio  does  not  equal  or  exceed  the  ratio  setting 
d,  which  is  a  variable  which  takes  on  a  value  greater  than  1.0,  then  declare 
an  erasure.  If  the  ratio  is  equal  to  or  greater  than  d,  then  declare  the 
symbol  corresponding  to  the  largest  correlator  output  as  the  transmitted 
symbol. 

The  ratio  rule  word  error  performance  for  analog  correlation  is  analyzed 
for  orthogonal  signalling  and  antipodal  signalling. 


For  the  present  analysis  of  the  ratio  rule  performance,  (4-2)  and  (4-3) 


hold  true. 


To  get  a  correct  symbol  decision,  r0  must  be  the  largest  correlator 

output,  and  it  must  exceed  the  next  largest  correlator  output  by  a  factor  d; 

r0  is  gaussianly  distributed  around  the  mean  energy  value  of  Es  with  variance 

E  N  /2.  The  other  correlator  outputs  are  similarly  gaussianly  distributed, 
s  o 

but  the  distributions  have  the  mean  value  of  (1/2)E  /N  .  the  probability  of  a 

so 

correct  decision  is 

i  ^oo  2  2 

Pc  e~(r_Es)  /2a  U-Qf(r/d-Es/2)/o]}31dr.  (4-22) 


to  error  occurs  if  any  of  the  31  correlators,  other  than  the  0’th 
correlator,  gives  rise  to  an  output  r . ,  i=l,2,...,31,  which  exceeds  the  second 
largest  correlator  output  by  a  factor  d.  the  probability  of  this  occurring  is 


» oo  2  2 

jsr-  [  e"  (r_Es/2)  /2a  {1-Q  [  (r/d-E  /2)  /a) } 30 {1-Q  [  (r/d-EJ  /a]  }dr . 
Y*it  o  J-*o  s  s  (4—23) 


To  find  the  probability  of  an  erasure  we  take  the  probability  of  a 
correct  decision  and  the  probability  of  an  error  and  express  the  probability 
of  an  erasure  as  follows: 


P 


E 


(4-24) 


Note  that  the  integration  of  a  gaussian  distribution  from  minus  infinity  to 
infinity  is  equal  to  one.  The  relationship  is  expressed  as 


(4-25) 


1 


-(r-E  )2/2o2 
e  s  dr. 


Using  (4-22),  (4-24),  and  (4-25),  we  get 

2  2 

PE  =  jTihfZ  e"(r_Es)  /2°  U-{l-Q((r/d-Es/2)/o)}31]dr  -  Pe.  (4-26) 

For  a  very  small  P£,  the  two-term  binomial  expansion  in  (4-13)  can  be  used  to 
approximate  (4-26)  as 

31  f  “  -lr- E 

PE  “  TSTo  L  e  s  Q(r/d-Eg/2)/a]dr  -  P#.  (4-27) 

P@  is  numerically  computed  using  (4-23),  and  Pg  is  numerically  computed 
using  (4-26)  for  high  P_  values  and  (4-27)  for  low  P_  values.  With  the 

b  E 

obtained  data,  P  vs.  p  for  constant  values  of  E  /N  are  plotted  in  figure 
4.5.  The  loci  of  constant  t  are  also  shown.  Superimposing  the  P^  curves  from 
figure  3.1  on  the  above  curves,  a  Pw  vs.  SNR  plot,  which  enables  us  to  predict 
the  detection  system  performance,  is  generated. 


Figure  4.5  shows  that  the  ratio  setting  of  d=1.06  yields  low  word  error 
probabilities,  we  find  that  the  ratio  setting  is  the  optimum  setting  for  all 
values  of  E  /N  shown  in  the  P  vs  P  graph.  If  d  is  varied  over  the  range 
1.03  <  d  <  1.10,  then  the  effect  on  the  word  error  probability  will  be 


minimal 


Figure  4.5.  Ratio  rule  orthogonal  signalling 
performance  curves. 


For  the  present  analysis  of  the  ratio  rule  performance,  (4-16)  and  (4-17) 
hold  true. 

The  probability  analysis  for  antipodal  signalling  is  essentially  the  same 
as  the  analysis  for  orthogonal  signalling,  but  the  statistical  mean  value  for 
r.  (for  i*l,2,...,31)  is  zero  for  antipodal  signalling.  Without  retracing  the 
derivations  we  will  simply  state  the  necessary  relationships  for  the 
performance  analysis  as  follows 

pc  =  [2  e'(r"Es}  2/2°2{l-Q[r/dal  }31dr  (4-28) 

pe  ■  e~r2/2°2 { 1-Q  (r/da]  }30{l-Q[(r/d-Es)l/o)  )dr  (4-29) 

PE  *  jTTSJ  _e~<r‘Es)  /2a  [l-U-0(r/da))31]dr-pe.  (4-30) 

Using  the  definition  in  (4-16),  (4-29)  and  (4-30)  can  be  numerically 
computed.  By  plotting  the  curves  as  was  done  for  orthogonal  signalling  we  are 
able  to  predict  the  probability  of  a  word  error. 

The  curves  derived  from  the  numerical  analysis  of  (4-29)  and  (4-30)  are 
shown  in  fig.  4.6.  We  see  that  the  ratio  setting  of  d*1.18  yields  low  word 
error  probabilities.  Hie  ratio  setting  is  optimum  for  all  values  of  Eg/No 
shown  in  the  Pp  vs.  P  graph.  If  d  is  varied  over  the  range  1.11  <  d  <  1.28, 
then  the  effect  on  the  word  error  probability  will  be  minimal. 


Figure  4.6.  Ratio  rule  antipodal  signalling 
performance  curves. 


4.3  Discussion  of  the  Results  of  the  Analog  Correlation  Analysis 


We  analyzed  the  performance  of  a  detection  system  using  the  threshold 
rule  and  of  a  detection  system  using  the  ratio  rule  for  orthogonal  signalling 
and  antipodal  signalling. 

Comparing  the  performance  results  of  the  threshold  rule  and  of  the  ratio 
rule,  we  find  that  the  optimum  ratio  setting  for  the  ratio  rule  is  relatively 
invariant  with  respect  to  changes  in  SNR,  whereas  we  see  the  optimum  threshold 
setting  for  the  threshold  rule  varying  with  changes  in  SNR.  When  the 
detection  system  performance  for  ratio  rule  is  compared  with  the  performance 
for  threshold  rule  the  ratio  rule  performance  is  slightly  better.  We  can 
conclude  that  for  the  case  of  analog  correlation,  the  ratio  rule  is  better 
than  the  threshold  rule. 

When  the  detection  system  performance  for  orthogonal  signalling  and  for 
antipodal  signalling  are  compared,  we  find  that  antipodal  signalling  achieves 
the  same  probability  of  a  word  error  as  orthogonal  signalling  for  6  dB  less 
than  the  SNR  required  for  orthogonal  signalling.  Therefore,  antipodal 
signalling  proves  to  be  better  than  orthogonal  signalling  in  terms  of 
detection  system  performance. 

In  this  chapter  we  analyzed  the  detection  system  performance  for  analog 
correlation.  As  seen  in  fig.  4.2,  the  symbol  detector  consists  of  32 
correlators,  each  holding  a  reference  signal  corresponding  to  a  32-ary  symbol. 
Each  correlator  produces  an  analog  output.  The  correlator  outputs  are  used  by 
the  decision  logic  device  to  make  a  symbol  decision. 
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CHAPTER  5 


DIGITAL  SIGNAL  DETECTION 


This  chapter  investigates  the  detection  system  performance  analysis  for 
the  case  of  digital  correlation.  For  this  case,  the  received  signals  are 
digitized  chi  a  chip-by-chip  basis. 

The  chief  advantage  to  the  use  of  digitization  is  that  it  simplifies  the 
hardware  design  and  implementation.  But  digitization  carries  the  penalty  of  a 
small  performance  degradation  when  compared  to  the  system  performance  for 
analog  signal  processing. 

A  block  diagram  of  the  digital  detection  system  is  shown  in  fig.  5.1. 

For  the  digital  detection  system,  a  chip  detector  precedes  the  bank  of  digital 
correlators.  There  are  32  digital  correlators,  each  corresponding  to  a  32-ary 
symbol. 

Each  symbol  is  represented  by  a  code  of  32  chips.  An  i'th  correlator 
correlates  the  code  from  the  chip  detector  with  its  reference  code  c.  and 
generates  the  quantized  correlator  output  rj,  where  r  ^-32, 30,..., 0,...,  30, 32 
and  i=0,l,...,31.  Upon  completion  of  the  digital  correlation  process,  the 
decision  logic  device  takes  the  symbol  correlator  outputs  and  decides  upon  a 
symbol  from  the  symbol  set 

U  *  Mi'  •••'  U31'  ensure}  (5-1) 


using  either  the  threshold  rule  or  the  ratio  rule.  The  symbol  decision  is 


sent  to  the  Reed-Solomon  decoder.  The  decoder  takes  the  31  symbol  decisions 
comprising  an  R-S  (31,15)  code  and  decodes  it  into  15  information  symbols. 


Figure  5.1.  Diagram  of  the  digital  detection  system. 


Hard  decisions  are  made  by  the  chip  detector  for  each  received  chip 
signal.  For  simplicity  of  analysis,  we  assume  that  each  hard  decision  yields 
a  value  of  1  or  -1,  depending  upon  the  received  chip  signal. 

Important  assumptions  are  made  about  the  chips  making  up  a  symbol.  The 
32  chips  representing  the  i'th  symbol  are  polar  (1,-1)  digits  comprising  a 
code  word  c . ,  where  i =0,1,..., 31.  Furthermore,  each  code  word  c.  has 
orthogonal  correlation  properties.  We  assume  that  when  a  code  word  is 
correlated  with  itself,  there  are  32  chip  agreements,  and  when  a  code  word  c^ 
is  correlated  with  c^  (for  i£j),  there  are  16  chip  agreements  and  16  chip 
disagreements.  Since  for  PSK  each  chip  agreement  yields  1  and  each  chip 
disagreement  yields  -1,  correlation  of  Cj  with  itself  will  yield  an  output  of 
32,  and  correlation  of  c^  with  c^  (for  ifj)  will  yield  an  output  of  0. 


The  decision  rules  developed  and  analyzed  for  analog  correlation  in  the 
last  chapter  will  be  analyzed  for  digital  correlation  in  this  chapter.  The 
results  of  the  digital  correlation  analysis  will  be  compared  with  the  results 
of  the  analog  correlation  analysis.  Because  of  the  quantization  of  signals  in 
digital  correlation,  we  expect  some  level  of  performance  degradation  when 
compared  to  the  system  performance  achieved  with  analog  correlation. 

For  this  chapter,  we  considered  antipodal  signalling  for  modulating  the 
chips.  For  orthogonal  signalling,  we  add  3  dB  to  each  constant  SNR  curve. 

5.1  Analysis  of  Digital  Signal  Detection 

Assuming  that  the  chips  comprising  a  code  word  are  transmitted  using 
antipodal  signalling  and  the  energy  in  the  received  signal  is  Es,  the 
probability  of  a  chip  error  occurring  at  the  chip  detector  is 

pce  “  o[VV32No]  (5-2> 

where  Eg/32  is  the  energy  of  each  PSK  modulated  chip  and  Nq/2  is  the  power 
spectra]  density  of  white  noise.  The  relationship  in  (5-2)  is  used  to 
formulate  various  probability  relationships  necessary  for  performance  analysis 
of  digital  signal  detection.  We  study  the  four  basic  probability  analyses 
fundamental  to  our  performance  analysis. 


Case  1;  The  probability  that  r^2r t  for  i^0. 

There  are  32  digital  correlators  in  the  digital  detection  system.  Since 
it  was  assumed  that  x0(t)  is  the  transmitted  signal  the  correlator  outputs  are 
seperated  into  two  cases:  r„  and  r.  (for  i=l,2,...,31).  Because  we  are 

fcJ  1 

dealing  with  orthogonal  codes,  for  noiseless  channel,  the  correlator  outputs 
have  the  values 

r0  =  32  (5-3a) 

rt  =  0  i=l,2,...,31.  (5-3b) 

We  will  refer  to  r,  for  which  i=l,2,...,31  as  r.,  and  r.  for  which  i=0  as  r_ 

X  11  v 

throughout  the  report. 

The  range  of  r.  depends  on  the  0  to  32  possible  chip  errors  in  a  32-chip 
sequence.  Figure  5.2  illustrates  all  possible  values  of  r^  given  k,  the 
number  of  chip  errors. 

We  explain  fig.  5.2  by  using  examples.  With  no  error,  r^  is  0  because  it 
is  the  sum  of  16  l's  (sign  agreements)  and  16  -l's  (si —  disagreements)  in  the 
chip-by-chip  correlation  of  c^  and  c^,  where  i#j.  With  one  error,  r^  can  be 
-2  when  it  is  the  sum  of  15  l's  and  17  -l's,  or  r.  can  be  2  when  it  is  the 
sum  of  17  l's  and  15  -l's.  We  can  do  the  same  calculations  for  other  k  values 
(the  number  of  chip  errors)  and  verify  the  results  in  fig.  5.2. 


The  probability  that  is  greater  than  or  equal  to  t  can  be  expressed  as 

b2-t/2j 

P(r.kt)  =  7  P[r,it)k]P(k)  (5-4) 

1  k=ft/2l  1 

where  (32-t/2j  is  the  floor  of  32-t/2f  or  the  greatest  integer  less  than  or 
equal  to  32-t/2,  and  ft/2~}  is  the  ceiling  of  t/2,  or  the  least  integer  greater 
than  or  equal  to  t/2.  The  probability  is  equivalent  to  the  sum  of  the 
probabilities  of  all  the  coordinate  points  located  at  or  above  the  line  r -=t. 


P(k)  is  the  probability  of  k  errors  occurring  in  a  32-chip  code  word. 

The  k  errors  in  a  chip  sequence  can  arrange  themselves  in  many  different  ways, 
requiring  a  combinatorial  analysis.  P(k)  is  expressed  as 

p<k>-  Cl)  p=e<1-t'ce>3Ml  (5-5» 

where  Pce  is  the  probability  of  a  chip  error  as  defined  in  (5-2),  and  (1-Pce) 
is  the  probability  of  an  error-free  chip.  Combining  (5-4)  and  (5-5),  we 
obtain 

|_32-t/2j 

P(r -it]  =  Z  P[r  .Jt|k] 

1  k* (t/2]  1 


Cl) 


P.Al-PJ 


32-k 


(5-6) 


1t»e  number  of  chip  errors  k  can  be  defined  as 


where  k^  is  the  number  of  chip  errors  which  cause  a  positive  contribution  to 
the  correlator  output  and  k.2  is  the  number  of  chip  errors  which  cause  a 
negative  contribution  to  the  correlator  output.  The  reader  is  referred  to 
fig.  5.2  to  note  that  chip  errors  cause  only  even  integral  contributions  to 
the  symbol  correlator  outputs.  Using  (5-7) ,  (5-6)  is  expanded  to  get 

32-k 


eiriit)  -  £  £  CjXk-iSp  P«<1H 

k*pt/2*|  k1=f(t+2k)/4]A 


P  )' 
ce 


L32t/2jl6  (16W  16\ 

Zr  _  L  r  W  W 

k=[t/2]  k  -f(t+2k)/4]1  1 


P  ( 1— P  ) 
ce '  ce' 


32-k 


(5-8) 


Equation  (5-8)  was  derived  by  summing  the  probabilities  of  the  coordinate 
points  occurring  at  or  above  the  horizontal  t  line  along  the  vertical  constant 
k  lines.  This  approach  required  us  to  break  the  expression  into  two  terms.  A 
more  simplified  equation  was  later  derived  when  the  probabilities  of  the 
points  occurring  at  or  above  the  horizontal  t  line  along  the  diagonal  constant 
k^  lines  were  summed.  The  expression  derived  is 


16 


Plrjjt] 


fjrt/2*j.i6.  .  16 

Z  .  _  I  'V  '>k‘k 

k1=ft/2]  k=kl  1 
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P  (1-P  ) 

ce  ce 


32-k 


(5-9) 


This  is  the  final  expression  describing  the  probability  that  r.  is  greater 
than  or  equal  to  t.  The  two  independent  variables  in  the  equation  are  t  and 


Case  2:  The  situation  in  which  r.2t 


We  would  like  to  find  the  probability  that  r0  is  greater  than  or  equal  to 
t.  Note  that  r0  is  32  when  there  exist  no  chip  errors.  But  with  each  chip 
error,  the  r0  value  is  decremented  by  2.  The  maximum  allowable  number  of  chip 
errors  in  a  32-chip  code  is  32,  in  which  case  the  r0  value  is  -32.  The 
probability  that  r0  is  greater  than  or  equal  to  t  is  written  as 


(16-C/2J 

I 

k=0 


P 

ce 


<1“Pce> 


32-k 


(5-10) 


The  two  independent  variables  in  (5-10)  are  t  and  Pce. 


Case  3;  The  case  in  which  r .=21. 

Referring  to  figure  5.2,  the  probability  that  r^  equals  an  even  integral 
value  21  is  the  sum  of  the  probabilities  of  the  points  occurring  coincident  to 
the  horizontal  line  r^«2I.  Along  such  a  line  k  and  kj  are  related  as 
4kj-2k*2I.  The  probability  is  expressed  as 


Case  4:  Ttie  case  m  which  r0=2I. 

Finding  the  probability  that  r0  equals  an  even  integral  value  21  is 
rather  straightforward.  Note  that  the  r0  value  in  the  absence  of  chip  errors 
is  32;  r0  and  k  are  related  as 

k  =  (32-r0)/2.  (5-1 2a) 

When  21  is  substituted  for  r_  in  (5-12a)  we  get 

0 

k  =16-1.  (5-1 2b) 


The  probability  that  r0  equals  21  is 


p(rr211  -  CD  <pc/<1-pc.> 


32-k 


(5-1 3a) 


But  vfien  16-1  is  substituted  for  k  we  get 


p(r<r2I>  -  (l«-3r)  <Pce,16'I(1-pce> 


16+1 


(5-1 3b) 


Having  obtained  the  expressions  for  Pfr^St],  P[r02:t],  P[r.=2I],  and 

P(r  =21]  we  can  proceed  to  evaluate  the  error  probabilities  and  erasure 
0 

probabilities  associated  with  the  threshold  rule  and  the  ratio  rule. 


5.2  Threshold  Rule  Analysis 

For  a  correct  symbol  decision  to  occur,  r0  must  be  greater  than  or  equal 
to  the  threshold  setting  t,  and  all  r j  (for  i=l,2,....,31)  must  be  less  than 
r0.  Since  the  correlator  outputs  are  even  integers,  they  can  be  represented 
as  21,  I  being  any  integer.  The  expression  for  the  probability  of  a  correct 
decision  is  given  by 

16 

Pr  -  Y  P(r  =2I)[l-P(r.2  2I)r\  (5-14) 

c  &  ft/2]  1 

For  an  error  to  occur,  one  of  the  r.  correlator  outputs  must  be  the 
largest  correlator  output  and  must  equal  or  exceed  the  threshold  t,  and  the 
rest  of  the  correlator  outputs,  including  r0,  must  be  less  than  the  largest 
correlator  output  value.  The  expression  for  the  probability  of  an  error  is 
given  by 

16 

Pe  -  31  £  P(ri=2I)[P(r.c2I)]30P(rc,*c2I).  (5-15a) 

I=ft/2]  1  1  B 

Putting  this  in  terms  of  the  previously  defined  terms,  we  obtain 
16 

P.  -  31  E  P(r.=2I)[l-P(r  »  2I)]30fl-P(r»2I)].  (5-15b) 

e  I- ft/21  1  1  0 

The  probability  of  an  erasure  is  formulated  by  simply  taking  the 


probability  of  a  correct  decision  and  the  probability  of  an  error,  which  are 
defined  in  (5-14)  and  (5-1 5b),  and  subtracting  them  from  1  as  follows 


PE  and  Pfi  are  numerically  computed  and  the  curves  are  plotted  to  predict 
the  digital  detection  system  performance  using  the  threshold  rule.  The 
performance  curves  found  in  fig.  S.3  show  that  the  threshold  should  be  set  at 
t*>12  to  yield  the  lowest  word  error  probability.  The  threshold  setting  can  be 
varied  in  the  range  of  10  <  t  <  14  without  substantially  affecting  the 
probability  of  a  word  error. 


For  a  correct  symbol  decision  to  occur,  r0  must  be  greater  than  or  equal 
to  d  times  the  largest  r.  value,  where  d  is  the  ratio  setting  which  can  take 
on  any  real  value  that  is  greater  than  1.  When  errors  are  introduced,  the  r0 
and  rj  values  must  be  even  integers  ranging  from  -32  to  32.  For  the  ratio 
rule,  we  approximate  our  analysis  by  assuming  that  the  probability  of  17  or 
more  chip  errors  occurring  in  a  code  of  32  chips  is  insignificant.  Therefore, 
we  restrict  our  analysis  to  the  case  where  only  16  or  fewer  chip  errors  occur. 
We  make  the  above  restriction  because  the  ratio  rule  performance  analysis 
becomes  complicated  when  negative  r0  values  are  introduced.  We  express  the 
probability  of  a  correct  decision  as 
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Pc  «  E0P(r0«2I)[l-P(ri2d(2I))]31.  (5-17) 

For  an  error  to  occur,  one  of  the  ri  correlator  outputs  must  be  the 
largest  correlator  output  which  equals  or  exceeds  a  value  that  is  d  times  the 
value  of  the  second  largest  correlator  output,  where  d  is  the  variable  ratio 
setting  which  takes  on  a  real  value  that  is  greater  than  1.  The  probability 
of  an  error  can  be  expressed  as 

16 

P  *  31  £  P(r.=2I)[l-P(r.2d(2I))]30[l-P(r.2d(2I)))  (5-18) 

e  1=0  1 

The  analysis  of  the  probability  of  an  erasure  is  not  easy  because  it 
involves  many  different  cases  of  an  erasure,  so  we  define  the  probability  of 
an  erasure  as  1  minus  the  sum  of  the  probability  of  a  correct  decision  and  the 


probability  of  an  error,  or 


1  -  P  -  P  . 
c  e 


(5-19) 


We  can  numerically  compute  the  P_  and  P  equations  and  plot  the  curves  to 

b  © 

predict  the  digital  detection  system  performance  vrtien  the  ratio  rule  is  used. 
The  performance  curves  are  shown  in  fig.  5.4.  The  optimum  value  of  d  is  1.25 
for  the  SNR  range  of  8  dB  to  11  dB.  It  is  not  critical  to  nMntain  this  value 
of  d  because  low  word  error  probabilities  can  be  obtained  for  a  broad  range  of 
values  of  d. 


1 


5.4  Discussion  of  the  Results  of  the  Digital  Correlation  Analysis 

Comparing  the  performance  results  of  the  two  decision  rules,  we  find  the 
optimun  ratio  setting  for  the  ratio  rule  and  the  optimum  threshold  setting  for 
the  threshold  rule  are  relatively  invariant  with  respect  to  changes  in  Es/Nq. 
This  was  not  the  case  for  the  threshold  rule  of  analog  correlation.  For  the 
digital  correlation  case  studied  in  this  chapter  we  find  no  distinct  advantage 
of  one  decision  rule  over  the  other. 

Recall  that  for  the  digital  correlation  case  the  chip  detector  outputs 
are  quantized.  This  means  that  the  correlation  values  are  quantized  into  even 
integral  values  ranging  from  -32  to  32.  The  quantization  of  the  correlator 
outputs  also  quantizes  the  values  that  the  threshold  setting  of  the  threshold 
rule  can  take  on.  The  threshold  setting  can  have  an  even  integral  value  which 
ranges  from  0  to  32.  However,  the  values  that  the  ratio  setting  for  the  ratio 
rule  can  have  due  to  quantization  are  not  easy  to  see,  but  we  do  find  discrete 
ratio  settings  which  do  affect  the  system  performance.  For  example,  the 
system  performance  is  not  affected  unless  the  ratio  setting  makes  a  transition 
through  any  of  the  following  values:  1.0,1.15,1.18,1.20,1.25,1.33,1.41.  Of 
course,  d  can  have  higher  values,  but  they  were  not  observed  because  they  fell 
out  of  the  range  of  d  values  of  our  interest. 

There  is  a  noticeable  performance  degradation  when  digital  correlation  is 
used  instead  of  the  analog  correlation.  In  terms  of  SNR,  we  see  less  than  1 
dB  performance  degradation  when  digital  correlation  is  used.  The  small 


degradation  in  the  system  performance  is  acceptable  when  we  consider  that  the 
design  and  implementation  of  the  digital  detector  is  simpler  than  the  analog 
detector. 


CHAPTER  6 


MONTE  CARLO  SIMULATION 


In  Chapter  5  the  detection  system  performance  for  various  cases  was 
analyzed.  The  derived  equations  were  numerically  analyzed,  and  the  data 
obtained  were  used  to  plot  the  performance  curves.  In  this  chapter  we  try  to 
simulate  the  operation  of  the  digital  detection  system  and  obtain  performance 
data  which  can  be  compared  with  the  data  derived  from  the  analysis.  The  Monte 
Carlo  simulation  method  is  used  to  process  signals  which  are  noise  corrupted. 
Antipodal  signalling  was  considered.  The  decision  rule  used  was  the  threshold 
rule. 


6.1  The  Simulation  Program 

The  digital  detection  system  shown  in  fig.  5.1  was  chosen  for  the  Monte 
Carlo  simulation.  The  operation  of  the  chip  detector,  the  bank  of  32 
correlators,  and  the  decision  logic  device  using  the  threshold  rule  are 
modeled  by  the  simulation  program. 

The  system  uses  32-chip  codes  which  have  orthogonal  correlation 
properties.  The  maximal-length  sequences  and  the  Welti  codes  are  investigated 
for  use  as  the  operating  codes.  The  performances  for  the  two  cases  are 
compared  to  each  other  and  to  the  performance  predicted  by  the  analysis  in 
Chapter  5.  We  will  later  see  that  the  simulation  using  the  modified 
maximal-length  sequences  and  the  simulation  using  the  Welti  codes  yield 


identical  performance  results,  but  the  simulation  results  do  not  sufficiently 
agree  with  the  results  obtained  from  the  analysis. 


The  simulation  repetitively  processes  x0  (t)  which  is  corrupted  with 
random  noise;  10,000  samples  of  32  detected  chips  are  collected,  and  the 
symbol  decisions  are  evaluated  to  find  the  relative  frequencies  of  symbol 
erasures  and  symbol  errors. 


6.1.1  Chip  Detector 

The  operation  of  the  chip  detector  is  modeled  using  the  random  number 
generator.  We  use  the  relationship 

Pce  ’  (6-1> 

to  relate  the  given  SNR  to  P  .  Hie  random  number  generator  is  used  to 
generate  independent  chip  errors  at  the  relative  frequency  that  approaches  the 
probability  of  a  chip  error  for  a  large  number  of  chip  samples.  For  10,000 
samples  of  a  32-chip  code,  it  is  typical  to  find  the  relative  frequency  of 
chip  errors  accurate  to  3  significant  digits  when  compared  to  the  given 
probability  of  a  chip  error. 


6.1.2  Digital  Correlators 

An  i'th  digital  correlator  takes  the  32  detected  chips  and  correlates  it 
with  the  reference  code  c^.  The  correlation  output  r.  can  take  on  an  even 
integral  value  ranging  from  -32  to  32.  As  stated  earlier,  x^  is  assumed  to  be 
transmitted  without  any  loss  of  generality.  Since  we  are  dealing  with 
orthogonal  codes,  under  the  noiseless  condition,  we  should  find  that 


r 

r 

r 


0 

1 

2 


<C0,C0>  ~  32 
<C0'°1>  =  0 
<C0,C2>  “  0 


r31  *  <c0,c31>  * 


0 


6.1.3  Decision  Logic  Device 

The  decision  logic  device  takes  the  32  correlation  outputs  and  makes  a 
symbol  decision  using  the  threshold  rule.  The  threshold  rule  makes  a  symbol 
decision  if  the  threshold  criterion  is  met,  but  if  the  threshold  criterion  is 
not  met  then  an  erasure  is  declared.  An  erasure  is  declared  when  the  received 
signal  has  unacceptably  degraded  in  signal  quality  due  to  the  noise  present  in 


the  transmission  medium 


The  threshold  rule  is  implemented  in  the  simulation  as  follows:  If  the 
largest  correlator  output  equals  or  exceeds  the  threshold,  then  declare  the 
syntool  corresponding  to  the  correlator  as  the  received  syntool.  If  either  the 
largest  correlator  output  does  not  exceed  the  threshold  or  if  we  have  a  tie, 
then  declare  an  erasure. 

Since  we  know  vtoich  syntool  is  sent  (we  assumed  that  is  sent  throughout 
the  report),  we  can  take  a  syiribol  decision  and  determine  if  the  decision  was 
correct  or  if  it  was  in  error.  By  counting  the  occurrences  of  errors  and 
erasures,  we  are  able  to  find  the  relative  frequencies  of  errors  and  erasures. 

6.2  Codes 

Codes  are  a  set  of  code  wards  used  to  generate  signals.  Bor  the  spread 
spectrum  ccmunication  system  that  we  examined,  each  32-ar/  syntool  is 
represented  by  a  code  word  of  length  32.  Although  the  actual  codes  used  are 
psuedo-orthogonal,  for  simplicity  of  analysis,  we  assume  codes  with  orthogonal 
correlation  properties. 

The  two  codes  investigated  are  the  maxima l-length  sequences  and  the  Welti 
codes.  The  codes  were  used  to  simulate  the  system  operation  and  obtain  the 
performance  data.  The  data  were  used  to  plot  the  performance  curves. 


6.2.1  Maximal-Length  Sequences 


In  this  section,  a  methodology  for  producing  maximal-length  sequences  of 
length  31,  each  having  the  periodic  correlation  function 


6{t) 


T=fl 

T*0 


(6-2) 


is  discussed.  The  maximal-length  sequences  are  then  modified  to  yield  codes 
of  32  chips  having  the  orthogonal  correlation  properties.  If  the  reader  is 
interested  in  a  general  discussion  of  binary  maximal-length  sequences,  then 
the  reader  is  referred  to  Bhargava,  Haccoun,  Matyas,  and  Nuspl  [7]. 

Maximal-length  sequences,  or  m-sequences,  can  be  generated  by  a  linear 
feedback  shift-register  generator  (LFSRG).  Given  an  m-stage  LFSRG,  a  set  of 
maximal-length  sequences  of  length  N  can  be  generated,  where 


N  =  2m  -  1. 


(6-3) 


Although  we  would  like  to  get  32-^ary  code  words,  (6-3)  tells  us  that  this  is 
not  possible.  Wte  can  take  the  next  best  choice  of  generating  31-ary  code 
words.  As  it  turns  out,  it  is  possible  to  modify  the  maximal-length 
sequences  so  that  we  can  obtain  32-ary  code  words  with  orthogonal  correlation 
properties. 

To  get  a  set  of  maximal-length  sequences  of  length  31,  a  5-stage  LFSRG 
must  be  used.  We  use  the  equation, 

h(x)  =  h0  +  hjX  +  h2x2  +  h3x3  +  h4x4  +  h5xJ  (6-4a) 
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where  h(x)  is  the  associated  polynomial  of  the  shift  register  with  feedback 


coefficients  (hg,h^,...,hj).  If  the  feedback  coefficients  are  defined  as 
hg=l,  hj«0,  hj^l,  h^=0,  h^=0,  h^=l,  corresponding  to  the  entries  in  the 
connection  vector  (101001)  for  a  5-stage  shift  register,  the  equation  can  be 
simplified  as 

h(x)=  1  +  x2  +  x5.  (6-4b) 

The  5-stage  LF5RG  modeled  by  (6-4b)  is  shown  in  figure  6.1.  The  initial  state 
is  assumed  to  be  (10000),  but  any  set  of  5  binary  digits  that  is  not  (00000) 
can  be  used.  The  output  from  the  last  stage  of  the  shift  register  is  used  to 
sequentially  extract  maximal-length  sequences. 


Figure  6.1.  5-Stage  Linear  Feedback  Shift-Register  Generator. 

A  sample  maximal-length  sequence  is  found  in  fig.  6.2a.  A  second 
maximal-length  sequence  is  found  in  fig.  6.2b.  Each  maximal-length  sequence 
generated  by  an  LFSRG  is  a  cyclical  shift  of  another.  In  order  to  have  a 
binary  (+1,-1)  chip  sequence,  we  convert  all  0's  into  -l's.  The  binary  (1,0) 
sequence  in  fig.  6.2a  becomes  the  binary  (+1,-1)  sequence  in  fig.  6.2c,  and 
the  binary  (1,0)  sequence  in  fig.  6.2b  becomes  the  binary  (+1,-1)  sequence  in 


0001010111011000111110011010010 

a.  A  sample  maximal-length  sequence  of  length  31. 

0010101110110001111100110100100 

b.  Another  maximal-length  sequence  of  length  31. 

-1-1-1  1-1  1-1  111-11  l-l-l-l  1111  1-1-1  1  1-1  1-1-1  1-1 

c.  A  binary  (+1,-1)  sequence  obtained  from  the  sequence  in  fig.  6.2a. 

-1-1  1-1  1-1  111-11  l-l-l-l  1111  1-1-1  1  1-1  1-1-1  1-1-1 

d.  A  binary  (+1,-1)  sequence  obtained  from  the  sequence  in  fig.  6.2b. 

1  l-l-l-l-l-l  1  1-1-1  1-1  11-1111  1-1  1-1  l-l-l-l  1-1-1  1 

e.  A  sequence  of  chip  correlation  outputs. 


Figure  6.2.  Maximal-Length  Sequences  of  Length  31. 


Thirty-one  maximal-length  sequences  are  shown  in  fig.  6.3a.  Each 
sequence  can  be  used  to  represent  a  31-ary  code  word.  To  generate  32-ary 
codes  with  orthogonal  correlation  properties,  we  append  the  digit  -1  to  each 
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of  the  31  sequences,  and  we  add  the  32nd  sequence,  which  consists  of  32  l's, 
to  the  set  of  31  sequences.  The  32-ary  orthogonal  codes  are  illustrated  in 
fig.  6.3b.  The  code  words  in  fig.  6.3b  have  orthogonal  correlation  properties 
as  follows: 

a.  When  a  code  word  c^  is  correlated  with  itself,  we  get  r.  *  32, 
vfriere  r4  is  the  correlation  output  of  the  i'th  correlator. 

b.  When  a  code  word  is  correlated  with  a  code  word  c^ 

(where  i^j),  we  get  r.  =  0. 


-1-1-1  1-1  1-1  111-11  l-l-l-l  1111  1-1-1  1  1-1  1-1-1  1-1 
-1-1  1-1  1-1  111-11  l-l-l-l  1111  1-1-1  1  1-1  1-1-1  1-1-1 


-l-l-l-l  1-1  1-1  111-11  l-l-l-l *1  1  1  1  1-1-1  1  1-1  1-1-1  1 
a.  Maximal-length  sequences  of  length  31. 


-1-1-1  1-1  1-1  111-11  l-l-l-l  1111  1-1-1  1  1-1  1-1-1  1-1-1 
-1-1  1-1  1-1  111-11  l-l-l-l  1111  1-1-1  1  1-1  1-1-1  l-l-l-l 


-l-l-l-l  1-1  1-1  111-11  l-l-l-l* 1  111  1-1-1  1  1-1  1-1-1  1-1 
11111111111111111111111111111111 

b.  32-ary  orthogonal  codes. 


Figure  6.3.  Two  Sample  Codes. 


The  codes  of  length  32  in  fig.  6.3b  are  used  in  the  Monte  Carlo 


simulation  to  produce  the  signals  represent inn  the  32-ary  symbols. 
Performance  data  are  obtained  for  10,000  samples  of  noise  corrupted  signals. 
The  performance  curves  derived  from  the  data  are  shown  in  fig.  6.4. 


Figure  6.5  shows  that  the  simulation  performance  is  worse  than  the  performance 
predicted  by  the  analysis  by  about  -1/2  dB.  Because  of  this  difference  in  the 
performance  results,  we've  investigated  the  use  of  another  orthogonal  code. 


6.2.2  Welti  Codes 

The  Welti  codes  18]  have  orthogonal  correlation  properties.  The  codes 
are  generated  as  follows:  Let  D*  represent  the  k'th  code  word  of  length  2l: 

d£  *  (  A  :  B  )  (6-6) 

where  the  word  is  divided  into  sets  of  lengths  2*  *.  From  this  word  the 
following  are  obtained: 


o£+1  *  (  A  :  B  :  A  :  B  ) 

(6-7) 

D^i  =  (  A  :  B  :  A  :  B  ) 

(6-8) 

where  x 

*  complement  of  x.  For  example, 

D0  *  <  1  « 
dJ  -  <  1  -1) 

Thus 

Dg  =  (  1  1  I’D 

D*  =  (  1  1-1  1) 

i. 


1  -1  1) 


and 


Dj  =  (  1  -1  1  1) 

r>3  =  (  l  -i  -l  -1) 

dJ  -  (  1  1  1  -l  1  1  -1  1) 

•  (  1  1  1  -1  -1  -1  1  -1) 

Djj  *  (  1  1-1  1  1  1  1-D 

Dg  *  {  1  1-1  1  -1  -1  -1  1) 

dJ  ■  (  1  -1  1  1  1-1-1  -1) 

Dg  »  (  1  — 1  1  1-1  1  1  1) 

»  (  1  -1  -1  -1  1-1  1  1) 

dJ*  (1  -1  -1  -1  -1  1  -1  -1) 


Similarly,  we  can  generate  words 

5  5 

which  generates  the  words  D0,  Dj,  ... 
simulation  program.  The  set  of  Welti 
6.6. 


of  length  32.  A  Welti  codes  generator 

5 

,  D32  is  implemented  in  the  Monte  Carlo 
codes  of  length  32  is  shown  in  figure 


Results  were  obtained  from  the  Monte  Carlo  simulation  using  the  Welti 
codes.  With  the  obtained  data,  Pp  vs.  P  for  constant  values  of  E_/N  were 
plotted  in  figure  6.7.  The  loci  of  the  constant  t  are  also  shown.  The  Pw 
curves  from  figure  3.1  were  superimposed  on  the  above  curves  to  generate  a  Pw 
vs.  SNR  plot,  which  enables  us  to  predict  the  detection  system  performance. 
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Figure  6.6.  Welti  codes  of  length  32 


The  simulation  results  for  the  Welti  codes  agreed  very  well  when  compared 
to  the  simulation  results  for  the  modified  maximal-length  sequences.  When  the 
simulation  results  are  compared  with  the  results  obtained  from  the  performance 
analysis,  we  find  that  the  simulation  results  deviated  by  an  average  of  0.5 


6.3  Discussion  of  the  Results  of  the  Simulations 

Monte  Carlo  simulation  was  used  to  verify  the  results  obtained  from  the 
digital  detection  system  analysis,  The  operation  of  a  digital  detection 
system  using  the  threshold  rule  for  its  decision  logic  device  was  simulated. 

sets  of  orthogonal  codes  were  investigated  for  possible  variations  in  the 
performance  results. 

When  the  performance  results  of  the  Monte  Carlo  simulation  using  the 
maximal-length  sequences  are  compared  with  the  performance  results  of  the 
Monte  Carlo  simulation  using  the  Welti  codes,  the  results  agreed  very  well  as 
expected  for  any  orthogonal  codes  of  length  32.  The  simulation  results, 
however,  were  found  to  deviate  from  the  results  of  the  analysis  by  an  average 
of  0.5  dB.  Further  investigation  revealed  that  the  deviation  of  the 
simulation  results  from  the  results  of  the  analysis  is  due  to  a  slight 
statistical  correlation  among  the  noise  corrupted  orthogonal  codes.  An 
in-depth  analysis  of  the  statistical  behavior  of  the  noise  corrupted 
orthogonal  codes  is  recommended  for  future  research  in  this  area. 


75 


CHAPTER  7 


CONCLUSION 


We  started  our  investigation  with  the  analog  detection  system.  There  we 
investigated  the  performance  of  the  system  which  uses  the  threshold  rule  and 
the  system  which  uses  the  ratio  rule.  “Hie  ratio  rule  was  found  to  have  an 
optimum  ratio  setting  which  was  invariant  with  respect  to  changes  in  the  SNR. 

We  next  investigated  the  digital  detection  system.  The  threshold  rule 
and  the  ratio  rule  were  found  to  have  the  parameter  setting  which  was 
invariant  with  respect  to  changes  in  the  SNR.  Compared  to  the  performance  of 
analog  detection,  the  system  performance  for  digital  detection  was  worse  by 
less  than  1  dB. 

Finally,  Monte  Carlo  simulation  of  the  digital  detection  system  using  the 
threshold  rule  was  conducted.  The  results  obtained  from  the  Monte  Carlo 
simulation  using  the  maximal-length  sequences  agreed  with  the  results  obtained 
from  the  Monte  Carlo  simulation  using  the  Welti  codes.  An  average  of  0.5  dB 
deviation  was  found  when  the  simulation  results  were  compared  with  the  results 
of  the  digital  detection  system  performance  analysis.  An  in-depth  analysis  of 
the  statistical  behavior  of  the  noise-corrupted  orthogonal  codes  is 


recommended  for  future  research  in  this  area. 
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LIST  OF  ACRONYMS  AND  SYMBOLS 


C  command  and  control 

command,  control,  and  communications 

e  number  of  symbol  errors 

E  number  of  symbol  erasures 

Ec  peak  chip-autocorrelation  value 

Eg  peak  symbol-autocorrelation  value 

e  .  minimum  value  of  e  given  S 

min 

e  maximum  value  of  e  given  S 

max 

Es/No  signal-to-noise  ratio 

FSK  frequency-shift  keying 

LFSRG  linear  feedback  shift-register  generator 

P  probability  of  a  correct  symbol 

v 

Pce  probability  of  a  chip  error 

Pe  probability  of  a  symbol  error 

P£  probability  of  a  symbol  erasure 

P(k)  probability  of  k  errors  occurring  in  a  32-chip 

code  word 

Pw  probability  of  a  word  error 

PSK  phase-shift  keying 

R-S  Reed-Solomon  (code  word) 

S  a  value  which  is  equal  to  2e  plus  E 

SNR  signal-to-noise  ratio 


